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Abstract 

The massive Schwinger model is studied, using a density matrix renormalisa- 
tion group approach to the staggered lattice Hamiltonian version of the model. 
Lattice sizes up to 256 sites are calculated, and the estimates in the contin- 
uum limit are almost two orders of magnitude more accurate than previous 
calculations. Coleman's picture of 'half-asymptotic' particles at background 
field 6 = 7r is confirmed. The predicted phase transition at finite fermion mass 
(m/g) is accurately located, and demonstrated to belong in the 2D Ising uni- 
versality class. 

PACS Indices: 12.20.-m, 11.15.Ha 
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I. INTRODUCTION 



The Schwinger model ]I|J2], or quantum electrodynamics in one space and one time dimen- 
sion, exhibits many analogies with QCD, including confinement, chiral symmetry breaking, 
charge shielding, and a topological vacuum f3)-|6|]. It is a common test bed for the trial 
of new techniques for the study of QCD: for instance, several authors [|7| P~0T have recently 
discussed new methods of treating lattice fermions using the Schwinger model as an example. 

Our purpose in this paper is twofold. First, we aim to explore the physics of this model 
when an external 'background' electric field is applied, as discussed long ago in a beautiful 
paper by Coleman ||. Secondly, we wish to demonstrate the application of density matrix 
renormalisation group methods fTTl , |T2"|l to a model of this sort, with long-range, non-local 
Coulomb interactions. 

Coleman || showed that the physics of the Schwinger model is periodic in 9 = 2nF/g, 
where F is the applied 'background' electric field, and g is the elementary charge. In the 
special case 9 = tt, some amusing phenomena occur. Whereas the "quarks" in this model 
are generally confined by a classical linear potential, at 9 = tt it is possible for single quarks 
to appear deconfined, i.e. move freely, provided they do not cross another quark ("half- 
asymptotic particles"). Coleman |J also demonstrated that for 9 = tt a phase transition 
must occur at some finite value of m/g, where m is the quark mass. These arguments are 
reviewed more fully in Section II. 

There have been very few attempts to verify these predictions numerically, that we 
are aware of. Hamer, Kogut, Crewther and Mazzolini |13| used finite-lattice techniques 
to address the problem. They calculated the ground-state energy and 'string tension' as 
functions of 9 and m/g. They located the phase transition at 9 = tt to lie at (m/g) c = 
0.325(20), with a correlation length index v = 0.9(1). They also attempted to estimate 
the chiral order parameter; but it was later pointed out WM that the chiral order parameter 



actually suffers from a logarithmic divergence at finite m/g in this model. Schiller and Ranft 
f23| used Monte Carlo techniques to locate the phase transition at (m/g) c = 0.31(1). 

Many different numerical methods have been applied to the Schwinger model in zero 
background field, including strong-coupling series [|T^-[TB|, finite-lattice calculations JTP|-|21 



Monte Carlo calculations fl22"|-f2~6f , discrete light-cone quantization [ 27j and the related light- 
front Tamm-Dancoff PE| , |2"9"|] and "fast-moving frame" [|30| techniques, a recently proposed 
"contractor renormalisation group" method |^TJ , and finally a coupled-cluster expansion [|32| . 



For a recent review, see Sriganesh et al. [21]. Analytic calculations of the mass spectrum 



have been carried out using mass perturbation theory [|^,|33}-|3"5| for small fermion masses, 
and weak-coupling expansions |6|, |36| , |2l| for large fermion masses. 

In this paper, we apply a new technique, namely the density matrix renormalisation 
group (DMRG), which has been used with great success [Tl],[l2| f° r lattice spin models and 
lattice electron models such as the Hubbard model. It was also recently applied to a simple 
one-particle potential model with asymptotic freedom by Martin-Delgado and Sierra |38] 



It might be questioned whether the DMRG can successfully handle a model involving long- 
range Coulomb interactions such as the Schwinger model. We perform some calculations 
in zero background field, to show that in fact the approach works extremely well. It gives 
accurate results for lattices of up to 256 sites, and provides estimates of the continuum limit 
which are around 50 times more accurate than previous calculations. 
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The layout of the paper is as follows. In Section II we review Coleman's discussion || of 
the Schwinger model in a background field, and the lattice formulation thereof. In Section III 
the DMRG technique is outlined. In Section IV our main results at background field 9 = it 
are presented. Mass gaps and order parameters are calculated, and the critical parameters 
at the phase transition are accurately estimated, showing that the transition lies in the 
universality class of the transverse Ising model in (1+1) dimensions. The 'half-asymptotic' 
particles behave in exactly the same way as 'kinks' or 'spinons' in the transverse Ising 
model, giving rise to some surprising effects. Finally, in Section V some results at 9 = 
are presented, and compared with previous results, both analytical and numerical. Our 
conclusions are summarized in Section VI. 



II. THE SCHWINGER MODEL IN A BACKGROUND FIELD 

A. Fermion formulation 

The Schwinger model Lagrangian density, in standard notation, is 

C = +${i$- g jL- m )$ (1) 

where ip is a 2-component spinor field, since there is no spin in one space dimension. The 
coupling g has dimensions of mass, so the theory is super-renormalisable. Using g as the 
scale of energy, the physical properties of the model are then functions of the dimensionless 
ratio m/g. The field strength term is 

= d^A v - d v A^. (2) 
The equations of motion are the Dirac equation 

(i - g 4-m)i/j = (3) 

and Maxwell's equations 

d fl F fU/ = gfyf1>- (4) 

Choosing a time-like axial gauge 

= (5) 

the field-strength tensor reduces to 

F 10 = —A 1 = E (6) 
where E is the 1-component electric field. Gauss' law becomes 

d x E = -d^ 1 = gf = # 7 V (7) 

The Hamiltonian becomes 
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(8) 



Gauss' law can be integrated to give 

E = gJ dxj°(x) + F, (9) 

showing that E is not an independent field, but can be determined in terms of the charge 
density j°(x), up to the constant of integration F, which corresponds to a "background 
field", as discussed by Coleman [El. We can think of the background field as created by 
condenser plates at either end of our one- dimensional universe. 

If \F\ > g/2, charged pairs will be produced, and separate to infinity, until the field is 
reduced within the range |F| < g/2, thus lowering the electrostatic energy per unit length 
(Fig. p]). Thus physics is periodic in F with period g, and it is convenient to define an angle 
9 by 

F 

9 = 2tt— . (10) 
9 

Then we can always choose 9 to lie in the interval [0, 2ir]. 

In the weak-coupling limit m/g — > oo, the vacuum contains no fermionic excitations, and 
the vacuum energy density eo is given purely by the electrostatic energy term (we ignore, or 
renormalise to zero, the energy of the Dirac sea). Hence 

_ ( \F 2 = g 2 9/ (8tt 2 ) (9<tt) 
° " 110? - F) 2 = g 2 (2n - 0) 2 /(8tt 2 ) (tt < 9 < 2tt). ^ 

Thus there is a discontinuity in the slope of the energy density, corresponding to a first-order 
phase transition, at 9 = n. In the strong-coupling limit m/g = 0, on the other hand, chiral 
invariance demands that the vacuum energy density remains constant as a function of 9 (see 
Sec. |IIB| ). Thus we expect a first-order transition at 9 = tt for large m/g, which terminates 



at a second order critical point at some finite (m/g) c , as illustrated in Fig. [| This behavior 



was demonstrated numerically some time ago by Hamer et al. [JL3], who located the critical 
point at (m/g) c = 0.325(20), with an associated critical index v = 0.9(1). 

Normally, charge is confined in the model: there is a 'string' of constant electric field (or 
flux) connecting any pair of opposite charges But Coleman points out that in the 

very special case 9 = it, or F = g/2, the peculiar phenomenon of "half-asymptotic" particles 
arises. In the weak-coupling limit, one can envisage the state shown in Fig. |3|. The electric 
field energy density is the same in between each pair of particles, and they can therefore 
move freely, as long as they maintain the same ordering, i.e. no pair of fermions interchanges 
positions. 



B. Boson formulation 



The one-dimensional fermionic theory can be mapped into an equivalent Bose form [p|,|37 
Some of the relevant mappings are 
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: ipip : < ► -cMN M cos(2^k4>) 

: z'V'75'0 : < — ► —cMN M sin(2A/7T0) 



n = n 



M 



f =: fyf$ : i > -^e^d^ (12) 
v vr 



(13) 



The Hamiltonian density in the charge zero sector can then be written 



^n 2 + -{d x (j)f + ^M 2 2 - cmM cos(2v^F0 - 9) 



Here A^m denotes normal ordering with respect to mass M, <p is the Bose field and IT is its 
conjugate momentum, and 
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M 2 = 9 —, (14) 

7T 



while 



c = k (15) 

and 7 = 0.5774 is Euler's constant. Note that 

M = (16) 

is Gauss' law, which identifies is proportional to the electric field. Finally, 9 is the back- 
ground field variable, as before. This form of the theory is most convenient for discussing 
the strong-coupling limit, m/g — ► 0. 

A nice discussion of the connection between 9 and chiral invariance is given by Creutz 
f5J. Consider the fermion mass term 

C m = mifjif). (17) 

If we now consider a chiral rotation 

i[) -> e idl5/2 i) (18) 



the mass term becomes 



where 



C m = m't/ji/j + m^iijj^ip (19) 



m' = m cos 9 

m5 = m sin 6 1 , (20) 

while the remaining terms of the Lagrangian density are naively invariant. In the Bose form 
of the theory, the transformed mass term maps into 

C m -> —otlMNm cos(2v^F(/) - 9) (21) 
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which is precisely the form ([13]), with the chiral rotation parameter 9 playing the role of the 
background field. 

It can be seen immediately from (|L|) that at m/g = 0, the Hamiltonian is independent 
of the background field variable 9 (i.e. is chiral invariant), and reduces to a theory of free, 
massive bosons, with mass M = g/\^ [1|J31- There is no sign of the "half-asymptotic" 
particles found in the weak-coupling limit at 9 = n. 

Coleman || gives a very neat semi-classical argument to show how the half-asymptotic 
particles arise in the Bose formulation. At 9 = ir, the Hamiltonian corresponds to an effective 
potential 

U{(j>) = ^M 2 <p 2 + cmMcos(2v^F0). (22) 

For m/g small, there is a unique vacuum at <p = 0. For m/g large, however, there are two 
vacua, located at = ±|-^/7r, and the symmetry <-> — suffers spontaneous breakdown. 
The two vacua correspond to background field values 9 = or 2n, or F = or g. Creutz |7j 
notes that in the broken symmetry phase the expectation value of (and hence sm(2y/7T(j))) 
should be non-zero, and therefore ipj^if) should make a suitable order parameter, by Eq. 
(0). 

Spontaneous breakdown of a Z2 symmetry implies that the critical point should belong 
to the universality class of the (1+1)D or 2D Ising model, with critical indices v = 1 and 



(3 = 1/8, which is consistent with the value for v found by Hamer et al. |L3 . 

For a single scalar field in two-dimensional space-time undergoing spontaneous symmetry 
breaking, there will exist time- independent finite-energy soliton solutions of the classical field 
equations passing monotonically from one ground state to the other. We may designate the 
increasing solution a "kink" and the decreasing solution an "antikink" (see Fig. [|). These 
are the half-asymptotic particles, in Bose language. Kinks and antikinks must alternate 
with each other when well separated, just as the half-asymptotic fermions do. 

At strong couplings, the mass of the single boson excitation at 9 = tt can be estimated 
semiclassically by fitting a harmonic oscillator to the effective potential at = 0. Hence one 
finds 



M 
9 



1 - v^e 7 - + O - 



= 0.564- 1.78 ^ +0^ , (23) 

where we have denoted this mass by A 2 , standing for the mass gap in the '2-particle', or 
electron-positron, sector. A crude linear extrapolation in m/g would give the boson mass 
vanishing at m/g = 0.317, quite close to the critical point found by Hamer et al. [p~3|. 



C. Lattice formulation 



We employ the Kogut-Susskind |39| , p!6| Hamiltonian spatial lattice formulation of the 
Schwinger model, with the fermions sited on a "staggered" spatial lattice. Let the lattice 
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spacing be a, and label the sites of the one- dimensional chain with an integer n. Define a 
single-component fermion field 4>(n) at each site n, and a link variable 

U{n, n + 1) = e ie{n) = e ' iagAl(n) (24) 
on each linkQ. Then the lattice Hamiltonian equivalent to (H) is 

N N 2 N 

H = -—Y / [^(n)e ie{n U(n + 1) - h.c] + m J2 (-l) n ^ (n)<f>(n) + — £ L 2 {n) (25) 

2° n=l n=l ^ n=l 

where the number of lattice sites N is even, and the correspondence between lattice and 
continuum fields is 

— 0(n) -> -A x (x) (27) 
<?L(n) - E{x). (28) 

The 7 matrices are represented by 

We use a "compact" formulation where the gauge field becomes an angular variable < 
6{n) < 2n on the lattice, and L(n) is the conjugate spin variable 

[9(n),L(m)]=i5 nm (30) 

so that L(n) has integer eigenvalues L(n) = 0,±1,±2, .... In the naive continuum limit 
a — > 0, the lattice Hamiltonian (E3) reduces to the continuum expression 
The Hamiltonian is transcribed to a dimensionless operator 

W = -^H = W + xV (31) 
ag 2 

where 

Wo = E lJ W+mE(- 1 )V(^W, (32) 

n n 

V = -iJ2 [<p f (n)e ie{n) 4)(n + 1) - h.c] (33) 



1 We trust the index (n) will differentiate the lattice fields <p(n) and 6(n) from the quite different 
fields and of Sec. |TT^ 
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In the lattice strong-coupling limit x — > 0, the unperturbed ground state |0) has 

L(n) = 0, 0t(n)0(n) = ^[l-(-l) n ], all n (35) 
whose energy we normalise to zero. The lattice version of Gauss' law is then taken as 

L{n) _ L{jl _ i) = 0t (n)0(n ) _ 1 [i _ (_!)»] (36) 

which means excitations on odd and even sites create =pl units of flux, corresponding to 
"electron" and "positron" excitations respectively. Eq. (^) determines the electric field 
L(n) entirely, up to an arbitrary additive constant a, which then represents the background 
field. Allowing a to be non-zero, the electrostatic energy term is modified to 

5> 2 W-E(^W+«) 2 - (37) 

n n 

The physics of the background field then matches precisely with the continuum discussion. 
If | a | > 1/2, charged pairs will be produced and separate to infinity to lower the electrostatic 
energy, and bring \a\ < 1/2. Physics is then periodic in a with period 1, and the background 
field variable is 



9 = 2vra. (38) 

In the weak-coupling limit m/g — > oo the vacuum contains no fermion excitations, and 
the ground state energy at x = is 

too = Na 2 (\a\ < 1/2) (39) 

corresponding to the 'zero-loop' state with {L(n) = 0, all n}; or 

uj = N(l-a) 2 (1/2<«<1) (40) 

corresponding to the 'one-loop' state with {L(n) = —1, all n}. 

At a = 1/2, these two states are degenerate; and half-asymptotic states consisting of 
alternating electrons and positrons exist, just as in the continuum model. Let us denote 
the two states ( p5Q and (|40| ) by ||) and | — |), according to the electric fields {L(n) + a = 
±|, all n} on the links. On a finite lattice at weak coupling, the eigenstates will be the 
symmetric and antisymmetric combinations 

|o> = ^[i|> + i-|>] (4i) 

and 

|0'> = ^[l§>-|-f>]- (42) 
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These states become degenerate in the bulk limit N — > oo, when spontaneous symmetry- 
breaking occurs. We shall denote the energy gap between them (or "loop gap") as A , 
standing for the O-particle sector. 



The lowest-energy single-particle excitation (Fig. |14]) is 



N 



|1> 



N ^ 



n = 1 
{n odd} 



N 

[n) J] e 

j=n 



I 21 



(43) 



while the lowest "2-particle" state is 



N 



N 

£ [0 f (n + l)e- 4e(n) 0WI|) + <t>\n + l)e* e(n) 0(™ + 2)| - \) 



(44) 



{n odd} 



We shall compute the energy gaps between these states and the ground state, Ax/g and 
A 2 /<7 respectively. 

We will also study two order parameters which can be used to characterize the phase 
transition at 9 = it. The first one is the average electric field 



1 

N 



(E(^)+a)>fl. 



(45) 



which in the weak-coupling limit m/g — > oo takes values ±1/2 for the zero-loop and one- 
loop states respectively. By the Feynman-Hellman theorem, it is proportional to the slope 
dE /da of the ground-state energy: a glance at Fig. || shows that at large m/g this slope 
undergoes a discontinuity at 9 = tt, while below (m/g) c it vanishes. 

The second order parameter is the axial fermion density suggested by Creutz 



(#7 5 ^/^) 
iJx 



N 



]T(-i) n f (n)0(n + 1) -h.c. ) 



where 



r 



z 7 y. 



(46) 
(47) 

(48) 



Now on a finite lattice there is no spontaneous symmetry breakdown, and the expectation 
values of the order parameters will remain identically zero. The remedy for this is well- 
known, however [40]. If the two lowest-energy states which become degenerate in the bulk 
limit iV — ► oo are denoted |0) and |0') respectively, then a simple 2x2 matrix calculation 
shows that the order parameter corresponding to an operator Q can be estimated as the 
overlap matrix element (0|Q|0') on the finite lattice. 

In zero background field a = 0, the 2-particle 'positronium' excited states of lowest 
energy in the lattice strong coupling limit are the 'vector' state |TS| 



1 

P= £ [<f)\n)e i9in) (f){n + 1) + h.c. |0> 

N n= l 



(49) 
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and the 'scalar' state 



1 

^) = 7^E iy(n)e*< B ty(n + 1) - h.c.l |0) 



(50) 



where |0) is the strong-coupling ground state (|35] 



D. Lattice spin formulation 



An equivalent lattice spin formulation can be obtained by a Jordan- Wigner transforma- 
tion fTBJ 



giving 



n 



n 



IT[^3( 



n[-^ 3 (/)]^ + H 



^o = E(^W + «) 2 + ^E(- 1 )^3H+ 



iV/x 



v = Y,[ 



a (n 



ie l9{n) a-(n + 1) + h.c. 



(51) 
(52) 



(53) 
(54) 



and 



L(n)-L(n-1) = -Mn) + (-!)"]■ 
The T 5 order parameter can be written in spin variables 



AT 



(E(-i: 



-+( n )e^ n )( 7 -(n+l)+h.c. 



This is the form which we used in the numerical calculations. 



III. METHODS 



A. The density matrix renormalisation group method 



(55) 



(56) 



Our results are based on the density matrix renormalisation group (DMRG) method 



11,12]. DMRG has been used primarily to study low dimensional quantum lattice systems 



in condensed matter physics, and is able to obtain with great accuracy quantities such as 
the ground and excited state energies, and correlation functions. One of the key features of 
DMRG is its ability to calculate these quantities for very large system sizes, particularly for 
systems with a low number of degrees of freedom per site. 

The method employed here is the "infinite system" DMRG method, as prescribed by 



White ||11|| , used both with open and periodic boundary conditions (OBC and PBC). Due 
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to the presence of the electric field on the links, and the differing nature of the odd and 
even-numbered sites, some modifications were made to the method, although nothing that 
changes the spirit of the DMRG. The quantities calculated are the ground state energies, 
mass gaps and order parameters. We use the form of the Hamiltonian given in (|53|) and 

Let us first look at how the presence of the electric fields affects the implementation of 
the code. For a particular spin configuration with OBC, if we specify the incoming electric 
field for the first site of the chain, then according to ( p5| ) the electric field for all the links can 
be deduced. This incoming field can be either L in = for the zero background field case, or 
L in = ±1/2 for the background field case. If PBC are imposed, then as there is no particular 
link to fix the electric field, we can have loops of electric flux extending throughout the ring. 
In the presence of a background field we simply add (or subtract) 1/2 from the values of the 
electric fields due to the spin configuration. A cutoff, or maximum loop value was chosen 
such that full convergence was reached to machine precision. A loop cutoff of [-5,5] was more 
than sufficient in most cases. 

A typical DMRG iteration is shown in Figs. [5] (OBC) and ^| (PBC). The chain (or ring) 
is split into two blocks and two sites, where blocks contain in general more than one site. In 
one DMRG iteration, we augment two sites to each block, so that the whole system grows 
by 4 sites each time. This is more convenient than the standard approach of adding one 
site at a time, due to the differing nature of the odd and even-numbered sites. For the PBC 
case this allows one to have the same type of site on each end of the block (i.e. both even or 
both odd), which considerably simplifies the bookkeeping of the electric fields on the links. 
In this case both blocks are identical, so that only one augmentation and density matrix is 
necessary. For the OBC case, we have a different situation. Again augmenting two sites at 
a time allows us to have the same type of site on each end of the blocks, but in this case 
the left hand block will always have an odd site on each end, while the right hand block will 
have an even site on each end. This requires two separate density matrices for each block, 
as they are not identical. 

We start with a 4 site superblock, and continue adding sites until we reach 256 sites. 
For PBC we keep a maximum of 125 states per spin sector per loop value in a block, which 
corresponds to a truncation parameter of approximately I ~ 1100, where I is the total 
number of basis states retained in a block. The truncation error for 256 sites can vary 
from 1 part in 10 7 to 10 13 depending on the number of excitations and the parameters. For 
OBC, again we keep a maximum of 125 states per spin sector, but as there are no electric 
field loops, here the truncation parameter is I ~ 300. The truncation error for 256 sites is 
typically 1 part in 10 16 . 

B. Accuracy of the DMRG calculations 

The accuracy of the DMRG calculation is strongly dependent on the parameter values x 
and m/g; in particular we expect the accuracy to worsen near critical regions. We find that 
OBC provide more accurate results than PBC, which is consistent with previous DMRG 
studies. We provide details of the accuracy for the worst case scenario, i.e. the parameters 
closest to the critical regions. This is not representative of the typical accuracies achievable 
with the DMRG, but gives us an upper bound on the errors. 
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We examine the point (m/g = 0.3; x = 100; 9 = tt; N = 256) with PBC, which represents 
the smallest lattice spacing used for this study, and lies in the critical region. Table | shows 
the convergence of the DMRG for these parameters as a function of the truncation parameter 
/, for the energies and order parameters. The ground state energy can be resolved to 1 part 
in 10 6 , while the two-particle gap is resolved to 1 part in 10 3 . The order parameters are 
subject to round-off errors because they involve the overlap of the two wavefunctions, and 
hence reduced accuracy is expected: here we have 3 figures. Fig. [7] shows the behaviour of 
the ground state energy density as a function of the density matrix truncation eigenvalue, 
which corresponds to the sum of the eigenvalues thrown away from the density matrix. We 
see a linear relationship between these quantities, which is confirmation that the DMRG is 
working correctly. Where possible, we obtain the final estimate by taking the results with 
the two largest truncation parameters, and performing a linear extrapolation to the y-axis. 
The error estimate is obtained by taking the difference in the extrapolant and the result 
with the largest /. 

Using OBC improves the accuracy of the DMRG calculation considerably. Table |I| shows 
the convergence of the same quantities with the same parameters as Table [I] , but with OBC. 
The ground state energy density is accurate to nearly machine precision, and the 2-particle 
mass gap to 1 part in 10 7 . Again due to round-off errors, we are limited in accuracy for the 
order parameters, so a similar accuracy to the PBC code is achieved here. A point against 
using OBC is that the finite-size corrections are much larger for this case. We use OBC to 
calculate the 1-particle gap and the order parameter estimates in the continuum limit, and 
PBC otherwise. 



IV. RESULTS AT BACKGROUND FIELD 9 = vr 

A. Analysis of critical behaviour 

We use finite-size scaling theory to estimate the position of the critical point, by cal- 
culating pseudo-critical points [^T| at each lattice size N, and each lattice spacing x. We 
demand that the scaled energy-gap ratio 

T> i i \ NA N (m/g) 

R " {mh) = (Jf-l)/W m / 9 ) (57) 

is equal to unity at the point m/g = (m/g)* N , which is the pseudo-critical point. Here 
A]\r(m/g) refers to the energy gap at some finite lattice size N. In practice we calculate 
R]sr{m/ g) for a cluster of five points straddling the pseudo-critical point, then use a polyno- 
mial interpolation to find (m/g)* N . The points were chosen at a spacing of A(m/g) = 0.02 
apart, the smallest reasonable spacing which ensures we cover the pseudo-critical point based 
on the work of Hamer et al. For this exercise we use the 'loop gap' A /g, which collapses 



to zero at the pseudo-critical point. 

Fig. |] shows three sample data sets for the pseudo-critical points calculated between 
x = 4 and x = 100. On a 1/N 3 plot these can be simply extrapolated with a quadratic, with 
errors in the vicinity of 1 part in 10 4 , for all x. Some numerical instability creeps in for the 
larger lattice sizes, however, where the change in the gap energy with lattice size N becomes 
so slow that round off errors become appreciable. The continuum limit is then estimated 
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from these bulk critical points, as shown in Fig. || A quadratic fit in 1/ y/x extracts the 
continuum limit a — > or x — > oo, which we estimate to be 



0.3335(2). (58) 



This is consistent with the previous estimate by Hamer et al. ITS] of (m/g) c = 0.325(20), or 



Schiller and Ranft [p3| , (m/g) c = 0.31(1), but with two orders of magnitude improvement 
in accuracy. 

Finite-size scaling theory |41[] also allows us to estimate the critical indices for the model. 



The theory tells us that the Callan-Symanzik "beta function" 

Mm/g) = A N (m/g) 

m/g (A N (m/g)-2(m/g)A' N (m/g)) 

scales at the critical point like ^((m/ 'g) c ) ~ N' 1 ^ as A" — > oo [Q. Hence we can extract 
the critical exponent v from the ratio 

iv(l- /~« m ^ ) . ), (60) 

since this should approach \/v as A" — > oo. Alternatively we can use the limiting behavior 

ln[p N ((m/ g y N )/p N -i((m/ g y N )] _1 

\n[N/(N - 1)] v 1 ' 

as A" — > oo. We call the first of these ratios flBTiD the 'linear' estimate, while the second is the 
'logarithmic' estimate. We calculate these ratios for all lattice sizes A" and lattice spacings 
x. Analysis of the data shows that the 'logarithmic' version is more stable numerically, and 
converges more quickly towards the bulk limit. We show an example of this convergence in 



Fig. 10 for a particular lattice spacing 1/y/x = ga = 0.45. We see monotonic convergence 
towards — 1/V — ► —1, until the same numerical errors begin to creep in as seen in Fig. [8]. 

A similar method was used to estimate the critical index (3. The order parameters 
and T% are expected to scale like ~ N~^/ u []42]| . 'Linear' and 'logarithmic' ratios are again 
constructed from these quantities, and analysed in the same way as for the "beta-functions" . 
Again the logarithmic ratios 

\n[T N ((m/gy N )/T N ^((m/gy N )] J_ 

In [Af/ (A" - 1)] v 1 ' 

seem to do better than the 'linear' ratios, in terms of numerical stability and convergence. 



Fig. [10] shows the estimate for (3/u, which uses the electric field order parameter. All final 



estimates for (3/v shown in Table |TXT| are calculated using = ((L + a))o- 

Table |TTT| displays our results for the critical exponents, for each different lattice spacing. 
We see essentially no variation in the exponents with lattice spacing, to within the accuracy 
of our calculations. Our best estimates for the critical exponents are thus 

v = 0.99(1) (63) 
P/u = 0.125(5). (64) 

These results provide reasonably conclusive evidence that the Schwinger model transition 
at 9 = 7i lies in the same universality class as the one dimensional transverse Ising model, 
or equivalently the 2D Ising model, with (u = 1, (3 = 1/8). 
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B. Behaviour in the continuum limit 



We now turn to estimating continuum limit values for the energy gaps and order pa- 



rameters. First of all, examination of the data points for 2-particle gaps in Fig. [TT] reveal 
that some extrapolation is necessary, especially for small lattice spacings, to obtain the bulk 
limit iV — > oo. In order to do so, we need to fit the data with an appropriate function of N. 
Now results from the last section revealed that the critical behavior of this model is closely 
related to that of the transverse Ising model. In that model the finite lattice values away 
from the critical point converge exponentially to the bulk limit, modulo half-integral powers 
of N [03 (see appendix). Hence we apply a form 




a + —= + — exp(-diV) + e (65) 



to the data, where a, b, c, d, e are fit parameters. Fig. [TT] shows that this form in fact fits the 
data very well, for all values of m/g. 

As a double-check of the results of the fit, we also apply a VBS sequence extrapolation 
routine H3| to the finite-lattice sequences, which provides an independent estimate of the 



bulk limit. Errors for the VBS estimates were obtained by examining the columns of the 
VBS extrapolants, and comparing the results with different VBS parameters a. Typically 
the VBS and the fit results agree to better than 1% accuracy, as long as we are not near the 
critical region where it becomes difficult to extract a reliable estimate. 

Fig. [12] shows the extrapolation to the continuum limit x — > oo, or a — > 0. We see that 
generally there is good agreement between the VBS extrapolations and the fit (|65|), except 
towards small lattice spacings where a discrepancy opens up. An extrapolation to the 
continuum limit is now performed by a simple polynomial fit in powers of \j \fx = ga. The 
double extrapolation therefore introduces quite large errors into our final results, compared 
to the original DMRG eigenvalues where the errors are at worst of order 0.1%. 



1. Loop energies 



Fig. [13] shows our results for the "loop" energy gap Ao/g for all values of m/g. We see 



that this gap vanishes at the same point as calculated in (p8[) , and is zero for all m/g > 
{m/g) a as predicted by Coleman H. The convergence to the bulk limit N — > oo for this 



case is much better than for the 2-particle gap illustrated in Fig. [TTJ, such that away from 
the critical region hardly any extrapolation is needed. This gives us good accuracy for the 
region m/g < 0.25, where the gaps are resolved to 4 figures. Near the critical region it 
becomes rather more difficult to extrapolate in a consistent manner. 



2. 2-particle gap 

We plot on the same figure results for the 2-particle gap A 2 /g. Again, the gap vanishes 
at the critical point (m/g) c , but on either side of the critical point there is a finite gap, 
and an almost linear behaviour with m/g. A linear fit through the points in the range 
m/g= [0, (m/g) c ] gives 
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A 2 (m\ 

— » 0.569 - 1.72 — , 66 
9 \9 J 

which agrees very well with the prediction made in (p3|). Note, however, that the behaviour 
is not exactly linear. 

The m/g = case is unique in the sense that we have a direct check on our results, since 
an analytic solution is known. For the massless case the background field has no effect on 
the physics of the model, since any background field is completely screened out. Hence all 
measurable quantities should be independent of the background field, and in particular the 
Schwinger boson mass |],|| 

— = w 0.56419. (67) 
9 

Our DMRG estimate gives A2/P = 0.57(1), which agrees with this result within errors. 



3. 1 -particle gap 

The 1-particle gap must be calculated using open boundary conditions, due to the mis- 
match in electric fields at either end of the chain. Fig. [14| shows the situation. A single 
charged particle in the middle will shift the electric field by ±g, and hence we have differing 
electric fields at either end. There is a further complication in the case with open boundaries, 
in that applying a background field of either a = 1/2 or —1/2 gives a different result for 
the ground state energy for a finite lattice. This is due to the 'staggered lattice' convention 
we have adopted, with 'electrons' appearing on odd sites, and 'positrons' appearing on even 
sites. The only consistent definition of the one-particle gap for a finite lattice is therefore 

A 1 /g = E 1 -±(E+ + E^) (68) 

where E\ is the 1-particle energy, Eq is the ground state energy with a = 1/2, and Eq that 
for a = -1/2. 

Apart from this complication the procedure in extrapolating to the continuum limit is 
the same as for the 2-particle state. Our results are shown in Fig. [13]. We see that the gap 
vanishes for m/g < (m/g) c , while for m/g > (m/g) c the 1-particle gap is very close to half 
the 2-particle gap. Once again, the behaviour is very nearly linear in m/g. 

The pattern of eigenvalues exhibited in Fig. |H| bears an extraordinary resemblance to 
that of the transverse Ising model [|P| , fl2]| , even down to the (almost) linear behaviour with 
m/g. In particular, we see that the energy of the 1-particle or 'kink' state vanishes at 
the critical point, and then remains degenerate with the ground state for (m/g) < (m/g) c . 
Assuming this degeneracy is exact, this indicates that a 'kink condensate' will form in the 
ground state for small mass, as discussed by Fradkin and Susskind f44|. It also indicates the 



existence of a 'dual symmetry' in the model, additional to that discussed in Section II, which 
is spontaneously broken in the low-mass region, and has not been explored hitherto. There 
should also be a 'dual order parameter' associated with this symmetry. By analogy with the 
Ising case, one would expect the 'dual' order parameter to be the expectation value of the 
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kink creation/destruction operator, and the 'dual' symmetry to correspond to inversion of 
this operator; but we have not verified this by explicit computation. 

The other notable feature of Fig. [13] is the degeneracy (within errors) between the zero- 
particle gap A /g and the 2-particle gap A 2 /g at small mass. Our physical pictures based 
on a weak-coupling representation do not give clear physical insight into this phenomenon. 

All the features just discussed are of course peculiar to the special case 9 = n. For 
instance, the 1-particle gap at any other value of 9 would be infinite. 



4- Order Parameters 

We may also obtain estimates for the order parameters T 5 and T a as functions of m/g 



Our results are displayed in Fig. |15| Both order parameters are zero, within errors, for 
m/g < (m/g) c . Near the critical region, particularly for m/g < (mjg) c it becomes quite 
difficult to obtain accurate estimates, but it appears that both order parameters turn over 
and drop abruptly to zero as the critical point is approached from above, consistent with 
the small exponent (3 = 1/8 found in the previous section. The axial density T 5 decreases 
steadily towards zero at large m/g, whereas T a approaches the expected value of 1/2. All 
our results for the gaps A /g, Ai/g, A 2 /g and both order parameters are shown in Table 
P"V| for future reference. 



V. RESULTS AT BACKGROUND FIELD 9 = 

The case of zero background field has been studied by many authors already, as outlined 
in the introduction. Our purpose here is to demonstrate how the DMRG can improve on the 
accuracy of existing results. The main quantities of interest are the "vector" and "scalar" 
state masses. The most accurate results to date are those of Sriganesh et al. [^TJ, who used 



numerical exact diagonalisation results together with a VBS extrapolation to obtain their 
final estimates. The largest lattice size calculated by these authors was N = 22, hence we 
expect to be able to do much better using our DMRG algorithm, which can go to much 
larger lattice sizes. 

For the 9 = case, we find that convergence with lattice size is much more rapid than 
for 9 = it, so much so that with N = 256, there is essentially no extrapolation necessary 



to obtain the bulk limit. Fig. [16| shows the data for m/g = 0, which is the exactly soluble 
case. For a particular lattice spacing x, we have very good convergence with lattice size, 
so that for x = 4 through to x = 100 we have 6 figure convergence. Even for the largest 
lattice spacing x = 400 used here, the gaps have converged to 4 figures. Our final continuum 
estimate is obtained by the method of linear, quadratic and cubic extrapolants in 1/ ^Jx~ as 



used in Ref. [21]. An example is shown in Fig. [17|. Our final estimate of m~ / g = 0.56419(4) 
for this case agrees extremely well with the analytic result flBTp. 

Table [V| summarises our results, which show between one and two orders of magnitude 
improvement in accuracy over the previous best estimates for small values of m/g. We must 
note, however, that for values of m/g > 1 we have little or no improvement in accuracy over 
previous results. To explain this, first we must note that the structure of the eigenvalue 
function in x shifts towards large x for large m/g. But at large x and large m/g there 
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are many "intruder" states below the vector state for finite N, artifacts of the finite lattice 
corresponding to states with no fermion excitations but with loops of electric flux winding 
around the entire lattice. This restricts the range of x that can be used, and hence the 
accuracy of the calculation. In a future calculation, a way of eliminating these intruder 
states must be found. 

A comparison of our results with previous works is shown in Fig. [Tj|. In the low mass 
region, we see excellent agreement between our results and the series expansion around 
m/g = 0. Our results are also fairly consistent with the results of Sriganesh et al, obtained 
by exact diagonalisation. The 'fast-moving frame' results of Kroger and Scheu seem to be 
consistently a little low in this region. In the large mass region the situation is reversed, 
as Kroger and Scheu obtain excellent agreement with the non-relativistic expansion, while 
both our results and those of Sriganesh et al. seem to be slightly high. We attribute his to 
the problem with intruder states discussed above. The Kroger and Scheu quasi-light-cone 
approach appears to remain more accurate in this weak-coupling region. 

Scalar state mass gaps m + /g may also be calculated using DMRG; here we merely 
demonstrate its applicability for one mass value m/g = 1.0, where we obtain m + / g — 2m/ g = 
1.118(1). This is again over an order of magnitude better that the previous best estimate of 
Sriganesh et al. who obtained for the same quantity 1.12(3). 



VI. CONCLUSIONS 

In this paper, we have demonstrated the application of the numerical DMRG technique of 
White and collaborators ITT] to a lattice gauge model, namely the massive Schwinger model. 
The results are in most cases nearly two orders of magnitude more accurate than previous 
calculations. Where exact diagonalisation can treat lattices up to 22 sites, DMRG gives very 
accurate results up to 256 sites. The long-range Coulomb interaction present in the model 
has proved no impediment to the DMRG technique: this is very likely connected with the 
fact that that the Coulomb interaction is screened fj\M, and the effective interactions are 
short-range. 

The most interesting results were obtained for the case of background field 6 = tt. We 
have performed a detailed study confirming the existence of the 'kinks' or 'half-asymptotic' 
particles predicted by Coleman |J, and have shown that there is a phase transition at 
(m/g) c = 0.3335(2) belonging in the universality class of the transverse or (1+1)D Ising 
model. The pattern of energy eigenvalues near the critical point bears a truly remarkable 
resemblance to the transverse Ising model, and points to some new physical effects, beyond 
those discussed by Coleman ||. 

In particular, for (m/g) < [m/g) c and 9 = tt, we find that the single-kink state is 



degenerate with the ground state. This points to the existence of a "kink condensate" |44 
in the ground state. It also points to the existence of a 'dual symmetry' and 'dual order 
parameter' in the model, which have not been discussed hitherto. We also find that the gap 
between the symmetric and antisymmetric 'loop state' combinations appears to be exactly 
degenerate with the 2-particle vector gap in this low-mass region. 

Some calculations were also carried out at zero background field, 6 = 0. The vector gap 
was obtained with great precision for small m/g; but at large m/g the accuracy was spoiled 
a little by finite-lattice "intruder" states. 
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The same methods should be applicable to other lattice gauge models in (1+1) dimen- 
sions. In higher dimensions, however, the DMRG technique does not have such a large 
comparative advantage over other techniques. Strenuous efforts are being made to develop 
improved DMRG algorithms for lattice models in (2+1) D; but nothing has even been at- 
tempted in (3+l)D, as far as we are aware. 
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VIII. APPENDIX 

Here we calculate analytically the form of the finite-size scaling corrections in the ID 



transverse Ising model. We follow the discussion by Hamer and Barber |42j| , hereafter 
referred to as I. 

The Hamiltonian is taken as 

M M 

H= X;(l-a 3 (m))-zX>i(m)<7i(m + l) (69) 

m=l m=l 



The model can be solved exactly by the methods of Schultz, Mattis and Lieb J45[. The 
mass gap on a finite lattice of M sites with periodic boundary conditions is 



where 



F(x, M) = 2(1 -x) + 2M[T 2M {x) - T M {x)} (70) 



T,(x)^ZM§) (71) 



and 



A(0) = [(l-x) 2 + 4xsm 2 9]V 2 . (72) 

For fixed x, as M — > oo, 

T M {x) = ^=T M (x) ~ c + (-)(! - * 2 ) 1/2 ^7T (73) 



where 



w c (x) = sinh \\—^) (74) 

Zy/ X 
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(Eq. (A1.16) of ref. I) and hence 

F(x, M) ~ 2(1 -x) -2tm 



'2(1 X 2 ) oMw 



xM 

which exhibits the leading finite-size corrections, as required. 
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TABLES 



TABLE I. DMRG estimates of the ground state energy density ojq/2Nx, the "2-particle" gap 
A2/3, and two order parameters, the mean field T a = ((L + a))o and axial density L 5 = (w/'75' i / , /<?}o 



as functions of I, the number of basis states retained per block. 


These results 


are for periodic 


boundary conditions, at x = 100, m/g = 0.3, 6 = ir, and ./V = 256 sites. 




I 


u /2Nx A 2 /g 


pa 


r 5 


244 


-.31676292 .26833 


.28068 


.29104 


324 


-.31677001 .25687 


.27598 


.28147 


416 


-.31677263 .25279 


.27417 


.27735 


555 


-.31677385 .25054 


.27319 


.27528 


742 


-.31677428 .24963 


.27273 


.27444 


932 


-.31677440 .24933 


.27266 


.27423 


TABLE II. As for Table |, with open boundary conditions. 


I 


uj /2Nx A 2 /g 


pa 


r 5 


107 


-.31611009314180 .19536745 


.256555 


-.257439 


162 


-.31611009382342 .19535871 


.256594 


-.257479 


199 


-.31611009385738 .19535765 


.256606 


-.257489 


244 


-.31611009386852 .19535733 


.256782 


-.257658 


327 


-.31611009387277 .19535723 


.256751 


-.257627 


396 


-.31611009387248 .19535720 


.256734 


-.257604 


TABLE III. 


Estimates for critical exponents 1/u and 0/v as 


functions of the lattice spacing 


parameter x = 


i/g 2 a 2 . 






X 


l/u 




P/u 


4.0 


0.99(1) 




0.126(5) 


4.938 


1.00(2) 




0.125(5) 


6.25 


0.99(2) 




0.125(5) 


8.163 


0.99(2) 




0.125(6) 


11.1 


1.00(4) 




0.126(6) 


16.0 


0.99(3) 




0.126(6) 


25.0 


0.99(3) 




0.127(6) 


44.4 


0.97(4) 




0.123(6) 


100.0 


1.0(1) 




0.12(1) 
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TABLE IV. Our results for the loop energy Ao/g, 1-particle gap Ai/g, 2-particle gap A2/g at 
background field 8 = tt. We also quote our results for the order parameters r a = {{L + a))o and 

r 5 = (^75^/5)0 • 



m/g 


A /g 


Ai/g 


A2/9 


pa 


r 5 


0.0 


0.5643(2) 


0.0± 10~ 6 


0.57(1) 


0.0 ± 10" 4 


0.0 ± 10" 4 


0.05 


0.4756(2) 




0.48(1) 






0.1 


0.3883(2) 


0.0 ± 10~ 6 


0.40(1) 


0.0 ± 10~ 3 


0.0 ± lO" 2 


0.15 


0.3020(5) 




0.30(2) 






0.2 


0.2173(5) 


0.0 ± 10~ 5 


0.23(4) 


0.00(2) 


0.000(5) 


0.25 


0.134(2) 




0.16(4) 






0.3 


0.05(2) 


0.0 ± 10~ 2 


0.03(7) 


0.0(3) 


0.0(2) 


0.4 


0.0 ± 10~ 3 


0.105(3) 


0.22(1) 


0.376(1) 


0.302(5) 


0.5 


0.0 ± 10" 4 


0.246(3) 


0.49(1) 


0.421(1) 


0.270(5) 


0.6 


0.0± 10~ 4 


0.3764(6) 


0.758(8) 


0.4430(5) 


0.238(5) 


0.7 


0.0 ± 10~ 4 


0.5020(2) 


1.006(4) 


0.4566(5) 


0.211(5) 


0.8 


0.0± 10~ 5 


0.6224(1) 


1.249(4) 


0.4657(5) 


0.189(5) 


1.0 


0.0 ± 10~ 5 


0.8530(4) 


1.711(4) 


0.4769(5) 


0.155(5) 



TABLE V. Comparison of bound-state energies for the "vector" state with previous works. The 
results of Sriganesh et al. and Crewther and Hamer were obtained through finite-lattice studies, 
while Eller et al., Mo and Perry, and Kroger and Scheu used light-cone or related methods to 
obtain their results. 



m/g 



This work 



Sriganesh 
et al. l2l 



Crewther & 
Hamer [09(1 



Eller 



Mo & 



Kroger & 






0.56419(4) 


0.563(1) 


0.56(1) 








0.125 


0.53950(7) 


0.543(2) 


0.54(1) 


0.58 


0.54 


0.528 


0.25 


0.51918(5) 


0.519(4) 


0.52(1) 


0.53 


0.52 


0.511 


0.5 


0.48747(2) 


0.485(3) 


0.50(1) 


0.49 


0.49 


0.489 


1 


0.4444(1) 


0.448(4) 


0.46(1) 


0.45 


0.45 


0.445 


2 


0.398(1) 


0.394(5) 


0.413(5) 


0.40 


0.40 


0.394 


4 


0.340(1) 


0.345(5) 


0.358(5) 


0.34 


0.34 


0.339 


8 


0.287(8) 


0.295(3) 


0.299(5) 


0.28 


0.29 


0.285 


16 


0.238(5) 


0.243(2) 


0.245(5) 


0.23 


0.24 


0.235 


32 


0.194(5) 


0.198(2) 


0.197(5) 


0.20 


0.20 


0.191 
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FIGURES 




FIG. 1. Creation of a charged pair alters the electric field by amount ±g in the intervening 
space. 




FIG. 2. Schematic plot of the vacuum energy density as a function of m/g and 9. The heavy 
line marks the first-order transition line where the energy density has a cusp, terminating at the 
second order critical point (m/g) c , where the discontinuity in slope with 6 goes to zero. 
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FIG. 4. "Kinks" and "Antikinks" are represented by transitions between the ground state 
solutions <p = ±ii/7r. Every transition corresponds to either an quark or an antiquark, according 
to the picture given in Fig. || 
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FIG. 5. A DMRG iteration using open boundary conditions. Two augmentations are needed 
to allow for the difference between odd and even sites. 
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FIG. 6. A DMRG iteration using periodic boundary conditions. 
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FIG. 7. Dependence of the estimated ground state energy density uiq/2Nx on the density matrix 
truncation eigenvalues, for x = 100, m/g = 0.3, 8 = it, using PBC. 
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FIG. 8. Pseudocritical points for various lattice spacings. 
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FIG. 9. Critical line in the m/g versus 1 / \[x = ga plane. Open circles are our present estimates, 
and squares are the previous results of Hamer et al. which are in good agreement. The dashed 
line is a quadratic fit to the data in ga. 
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FIG. 10. 'Logarithmic' ratio estimates of critical indices — 1/v and —(3/v for lattice spacing 
l/\fx = ga = 0.45. Quadratic fits in 1/N provide the bulk extrapolations. We estimate here 
1/v = 1.00(2) and 0/v = 0.125(5). 
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FIG. 11. Bulk extrapolations for the 2-particle gap f^ild for Tn/g = 0.0 and 9 = ir. Dashed 
lines are the fits to the data, according to (pq) 
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FIG. 12. Continuum extrapolations for the 2-particle gap for m/g = 0.0 and 9 = tt. Data sets 
obtained through separate VBS and fit extrapolations. 
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FIG. 13. Final estimates for gaps in the 0-particle, 1-particle and 2-particle sectors at 6 = n. 
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FIG. 14. The electric fields at either end of the chain in the presence of a single particle are 
g/2, —g/2 respectively. This precludes periodic boundary conditions due to the mismatch of electric 
fields. 
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FIG. 15. Order parameters T a = ((L + a))o, T 5 = (i^5ip/g)o near the critical region. Dashed 
lines are merely to guide the eye. 
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FIG. 16. "Vector" mass gaps m~ jg for m/g = 0, finite lattices iV = 10-256, and various lattice 
spacings. 
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FIG. 17. An example of our extrapolation procedure for the "vector" mass gaps at m/g = 0. 
Circles, squares and triangles show linear, quadratic and cubic extrapolants respectively. Dashed 
lines show the upper and lower bounds for our final estimate. Here we estimate m~ /g = 0.56419(4). 
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FIG. 18. Comparison of our results for the "vector" state binding energies E\/g with other 
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Scheu |3(| . The results of Vary, Fields and Pirner |?3| and Adam were used for the expansion 
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Sriganesh et al. [21]. 



38 



